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H. ABSTRACT 

SUMMARY 

The  electron  density  distribution  in  the  shock  layer  of  re¬ 
entry  bodies  in  hypersonic  flight  is  determined  by  many  factors. 
Among  these,  the  effect  of  a  nonzero  angle-of-attack  of  the  body 
to  the  free  stream  is  an  important  one.  This  report  considers  the 
three  dimensional  boundary  layer  oYer  a  sphere-cone  body  at  hyper¬ 
sonic  fligiit  and  a  nonzero  angle-of-attack  at  conditions  where  the 
air  is  a  mixture  of  dissociated  and  ionized  components  reacting 
with  each  other  at  finite  chemical  rates.  The  equations  of  motion 
for  the  three  dimensional  flow,  are  formulated  and  method  of  attack 
by  an  implicit  finite  difference  scheme  is  utilized  for  their  solu- 
tiohs.  Edge  conditions  with  partial  swallowing  for  the  boundary 
layer  are  obtained  directly  from  the  inviscid  equations  of  motion 
which  are  solved  by  an  iterative  finite  difference  predictor- 
corrector  method.  Finally,  illustrative  numerical  results  of  the 
three  dimensional  flow  effects  are  presented  for  flow  over  the 
AFCRL  Trailblazer  II  Body. 
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SUMMARY 


The  electron  density  distribution  in  the  shock  layer  of  re-entry  bodies 
in  hypersonic  flight  is  determined  by  many  factors.  Among  these,  the  effect 
of  a  nonzero  angle-of-attack  of  the  body  to  the  free  stream  is  an  important 
one.  This  report  considers  the  three  dimensional  boundary  layer  over  a  sphere- 
cone  body  at  hypersonic  flight  and  a  nonzero  angle-of-attack  at  conditions 
where  the  air  is  a  mixture  of  dissociated  and  ionized  components  reacting  with 
each  other  at  finite  chemical  rates.  The  equations  of  motion  for  the  three 
dimensional  flow  are  formulated  and  a  method  of  attack  by  an  implicit  finite 
difference  scheme  is  utilized  for  their  solutions.  Edge  conditions  with 
partial  swallowing  for  the  boundary  layer  are  obtained  directly  from  the  in- 
viscid  equations  of  motion  which  are  solved  by  an  iterative  finite  difference 
predictor-corrector  method.  Finally,  illustrative  numerical  results  of  the 
three  dimensional-  flow  effects  are  presented  for  flow  over  the  AFCRL  Trailblazer 
II  body. 
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I.  Introduction 


The  production  of  electrons  at  hypersonic  speeds  in  the  shock  layer 
of  a  blunt  body  is  a  direct  function  of  the  conversion  of  kinetic  energy 
of  flight  to  heat  which  is  dependent  principally  on  the  aerodynamic  flight 
configuration.  At  the  high  temperatures  created  electrons  are  obtained 
from  the  various  energy  exchanges  in  the  hot  dissociated  air  about  the 
body.  The  existence  of  a  layer  of  electrons  about  re-entry  body  leads 
to  several  detrimental  effects;  one  of  the  principal  effect  is  the  attenua¬ 
tion  of  electromagnetic  signals  propagating  through  the  layer  leading  to 
telemetry  blackout.  A  number  of  studies  has  been  made  of  this  problem. 
These  studies  have  involved  both  flight  measurements  of  the  severity  of 

the  problem  (2)  an(j  theoretical  predictions  of  the  degree  of  ioniza- 

(3),  (4),  (5) 

tion  in  the  shock  layer  with  comparison  of  flight  data  with 

(1)  (2) 

theoretical  results  .  Good  agreement  was  found  between  theoretical 

prediction  and  flight  measurement  of  electron  density  about  the  Trailblazer 
sphere-cone  body  at  hypersonic  speeds  at  applicable  conditions. 

The  prediction  of  electron  density  about  hypersonic  re-entry  vehicles 
in  the  viscous  region  near  the  body  was  given  in  reference  (3)  and  (4)  for 
the  sphere  cone  body.  Numerical  results  for  the  AFCRL  Trailblazer  II 
vehicle  at  hypersonic  speed  were  given.  These  results  were  made  princi¬ 
pally  for  bodies  at  zero  angle-of-attack.  Since  the  gaseous  environment 
is  composed  of  hot  dissociated  species  internal  energy  transfers  become 
important.  The  effect  of  energy  distribution  in  the  various  energy  modes, 
in  particular,  the  vibrational  relaxation  process  may  change  the  production 
of  electrons.  A  study  of  these  effect  and  surface  catalyticity  on  the 
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electron  density  distributions  was  considered  in  reference  (5)  with 

numerical  results  for  the  AFCRL  Trailblazer  vehicle. 

Flight  of  the  Trailblazer  vehicle  has  shown  that  the  angle-of- 

attack  is  of  the  magnitude  of  10  to  15  with  respect  to  the  free 

stream.  The  asymmetry  of  the  flow  introduced  by  this  flight  angle 

lead  to  significant  changes  in  the  flow  field  characteristics  of 

velocity,  temperature,  and  concentration  of  the  dissociated  and  ionized 

air  species.  The  magnitude  of  the  effect  of  the  non-zero  angle-of- 

attack  is  usually  determined  by  the  relative  sizes  of  the  angle-of- 

attack  and  the  half-cone  angle.  If  oC  is  the  angle-of-attac k  and 

0  the  cone  half-angle  then  small  angle-of-attack  effects  are  re- 
c 

stricted  to  values  of  (  oC  /  0)^1.  For  the  Trailblazer  vehicle  the 

c 

cone  angle  is  9°  so  that  o4  / ©c  >  1;  that  is,  the  effect  of  angle-of- 
attack  is  large.  Observation  of  some  of  the  effects  of  large  angle-of- 
attack  has  been  made.  The  pattern  for  flow  over  a  blunt  cone  at  non¬ 
zero  angle-of-attack  has  been  obtained  from  experiments.  Laboratory 

(6) 

experiments  by  Tracy  have  shown  that  the  viscous  layer  on  the 

leeward  side  of  a  circular  cone  grows  exponentially  to  the  point  of  flow 

separation  and  subsequently  alters  the  inviscid  flow.  Later  experiments 

(7) 

by  Stetson  in  1971  gave  additional  details.  Stetson's  measurement 
indicate  that  there  are  formed  symmetrical  helical  vortices  with  no  re¬ 
verse  flow. 

In  general,  the  non-zero  angle-of-attack  generates  two  effects  which 
are  important  to  the  distribution  of  electron  density  over  a  blunt  coni¬ 
cal  body.  First,  the  pressure  is  increased  on  the  windward  side.  This 
increase  can  be  large.  Examples  of  this  increase  can  be  seen  from  Stetson's 
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(7) 


experimental  measurement  of  the  flow  over  a  5.6  half-angle  blunt 

cone.  For  a  bluntness  ratio  of  0.1  and  for  X/R^  =  5  (X  is  the  wetted 
length  measured  from  the  tip,  Rjq  is  the  nose  radius),  the  surface  pres¬ 
sure  ratio  on  the  most  windward  meridian  plane  has  been  increased  by 
a  factor  of  2  as  the  angle-of-attack  is  increased  from  zero  to  10°. 

The  maximum  Trailblazer  vehicle  X/R^;  is  4.25  and  its  bluntness  ratio 
is  0.664  and  its  cone  half-angle  is  9°.  Thus  a  factor  of  2  to  3  in¬ 
crease  in  pressure  would  be  expected  on  the  windward  meridian  plane. 

This  surface  pressure  distribution  decreases  in  value  as  the  flow 
sweeps  toward  the  leeward  plane;  at  this  latter  plane  the  pressure  drops 
below  the  values  at  zero  angle-of-attack.  Thus  the  increase  of  pressure 
will  lead  to  a  noticeable  effect  on  the  electron  density  by  the  increase 
in  production  rate  of  electrons. 

A  second  effect  is  the  divergence  of  surface  streamlines  from  the 

(8) 

inviscid  streamlines.  Rakich,  et  al  has  made  an  oil-film  experi¬ 
ment  of  the  surface  streamlines  showing  the  turning  of  the  viscous 
streamlines  near  the  body  relative  to  the  inviscid  streamlines  as  shown 
on  Figure  1.1.  The  significance  of  the  cross  flow  can  be  appreciated  by 
reference  to  Figure  1.1  where  the  increasing  divergence  in  the  directions 
of  the  boundary  layer  and  inviscid  flow  streamlines  downstream  is  shown. 
The  region  of  influence  as  noted  in  the  Figure  outlines  the  entire  region 
influenced  by  the  original  planar  meridian  streamlines.  The  curvature 
of  the  fluid  streamlines  near  the  wall  becomes  much  different  from  that 
of  the  inviscid  streamlines  due  to  the  cross  flow  generated  by  a  non¬ 
zero  angle-of-attack.  Species  concentrations  are  changed  locally  by 
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convection  following  the  increased  divergency  of  streamlines. 

The  mathematical  characteristics  of  the  three  dimensional  bound¬ 
ary  layer  still  remain  parabolic.  However,  a  large  difference  appears 
with  the  addition  of  the  one  coordinate.  Characteristics  such  as 
momentum,  energy,  and  species  concentrations  are  diffused  normal  to 
the  surface  and  convected  by  streamlines  downstream  into  the  zone  of 
influence  shown  on  Figure  1.  1.  The  electron  concentration  generated 
at  the  stagnation  region  spread  over  the  windward  side  and  then  to  lee¬ 
ward  side  of  the  cone.  Thus  the  second  effect  of  non-zero  angle-of- 
attack  is  the  divergency  of  the  surface  streamlines  toward  the  lee¬ 
ward  plane.  This  results  in  convection  of  electrons  away  from  the  wind¬ 
ward  planes.  Consequently,  one  can  say  that  the  net  effect  of  a  non¬ 
zero  angle-of-attack  is  to  generate  more  electrons  and  to  transport 

•  v 

them  out  of  the  meridian  plane  to  all  other  parts  of  the  flow  field. 

This  report  is  concerned  with  the  continuation  of  the  previous 

(3,4,5) 

studies  to  determine  the  electron  density  distribution  over  blunt 

sphere  cone  bodies  at  hypersonic  speed  to  non-zero  angle-of-attack  of  the 
body  to  the  free  stream.  The  behavior  of  the  three  dimensional  laminar 
nonequilibrium  boundary  layer  on  blunt  conical  bodies  is  considered  with  em¬ 
phasis  on  chemical  nonequilibrium  aspect  of  the  air  species  and  in  part¬ 
icular,  the  production  of  electron.  As  in  the  previous  studies  diffusion 
of  the  dissociated  and  ionized  air  species  with  finite  rate  chemistry  is 
also  emphasized.  Numerical  results  are  given  for  the  AFCRL  Trailblazer 
body. 
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In  Section  II  the  equations  of  motion  are  given  for  the  flow  of  a  three 
dimensional  boundary  layer  in  an  orthogonal  curvilinear  coordinate  system 
appropriate  for  a  sphere-cone  body  at  angle-of-attack.  These  equations  are 
valid  for  a  multi-component  mixture  of  dissociated  and  ionized  gas  species 
which  may  be  reacting  with  each  other  with  finite  chemical  rates.  Viscosity, 
heat  conduction,  and  diffusion  are  the  principal  transport  phenomena  included 
in  the  formulation.  The  boundary  layer  developing  over  the  body  is  acted  upon 
by  the  outer  inviscid  flow  which  provides  the  edge  conditions  under  which  the 
boundary  layer  develops.  One  complication  as  considered  in  reference  (3)  is 
the  swallowing  of  the  streamlines  with  different  entropy  downstream  of  the 
stagnation  region. 

Section  III  contains  the  inviscid  flow  analysis;  considerations  are  given 
to  formulating  the  appropriate  set  for  the  three  dimensional  boundary  layer  edge 
conditions.  A  numerical  method  of  solution  is  provided  for  the  three  dimensional 
flow  of  a  chemical  reacting  gas.  Numerical  results  for  the  edge  conditions  are 
discussed  in  Section  V  and  are  given  for  a  air  mixture  with  finite  reacting  chem¬ 
istry.  Section  IV  summarizes  the  chemical  rate  constants,  thermodynamic,  and 
transport  properties  for  a  dissociating  and  slightly  ionized  air  mixture  at  high 
temperature  which  are  used  in  this  study. 

This  report  contains  the  equations  utilized  in  the  study  and  numerical  re¬ 
sults  for  boundary  layer  viscous  flow  at  the  stagnation  line  for  bodies  at  non¬ 
zero  angle-of-attack.  Numerical  results  are  given  for  the  edge  conditions  due 
to  the  three  dimensional  inviscid  flow.  Further  numerical  results  will  be  given 
in  a  subsequent  report. 
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II.  Equations  of  Motion  for  Three  Dimensional  Laminar  Boundary  Layer 

The  problem  of  concern  is  the  flow  of  a  three  dimensional 

boundary  layer  over  a  body  at  angle-of-attack  at  hypersonic  speeds 

with  transport  processes  and  chemical  reaction  between  dissociated 

and  ionized  air  species.  Since  species  of  different  molecular  weight 

are  present  diffusion  and  convection  of  each  species  dominates  the 

electron  density  at  any  point  of  the  body.  To  date  all  studies  of 

three  dimensional  boundary  layer  theory  have  emphasized  non-reacting 

perfect  gases.  These  have  been  reviewed,  for  example,  in  the  article 

(9)  (10) 

by  Mager  and  a  more  recent  one  by  Eichelbrenner 

It  is  convenient  to  use  a  set  of  orthogonal  curvilinear  coordinates 
where  the  surface  of  the  body  is  normal  to  one  of  the  coordinates.  With 
this  normal  coordinate  the  concept  of  a  thin  viscous  layer  developing 
over  the  surface  can  be  utilized  for  obtaining  the  boundary  layer  equa¬ 
tions  from  the  Navier-Stokes  equations  of  flow.  With  this  set  of  coordi¬ 
nates  the  line  element  assumes  the  form. 

dsx  =  + £5 d z*'  (2-1) 

where  x  and  z  are  the  surface  coordinates  and  y  is  normal  to  the  body 
surface  and  the  metric  coefficients  h^,  h2>'  and  I13  are  functions  only 
of  x  and  z  when  the  boundary  layer  is  thin  compared  to  the  principal 
radii  of  curvature  of  the  surface.  By  setting  h2  =  1,  y  becomes  the 
normal  distance  to  the  surface  and  it  remains  to  choose  the  orthogonal 
surface  coordinates.  It  is  convenient  for  this  study  to 
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(11) 

utilize  the  system  of  geodesic  coordinates  (Moore  )  which 

simplifies  the  scale  factor  further  to  hj_  =  1  and  the  second  one 

an  arbitrary  scale  function  h^  =  r  (x,  y).  For  rotationally 

symmetric  bodies  such  as  sphere  cones  x  is  measured  along  the  body 

generator,  z  the  meridional  angle,  and  h^  =  r  (x)  is  the  radius 

of  the  body.  Thus  the  only  non-unit  metric  coefficient  is  r  (x) . 

This  system  of  coordinate  is  the  same  as  that  used  in  the  zero  angle- 
(3) 

of-attack  study  .  For  any  other  more  general  body  the  set  of 
orthogonal  curvilinear  coordinates  can  be  obtained  by  a  selection  of 
one  family  and  then  obtaining  the  orthogonal  projections  by  numerical 
computation. 


In  this  set  of  coordinates  where  x  is  along  the  meridian  curve, 
z  is  the  meridional  angle,  and  y  is  the  distance  normal  to  the  surface, 
The  equations  for  the  conservation  of  momentum,  energy,  and  species 
for  the  boundary  layer  are  respectively  given  by: 


p(dU  „  w\ '  ) _ ,  d  fndU) 

J  L  dt  7€  J  ~  7k  +  (r 


(2.2) 


(2.3) 


d  t 


~d7t 


1st 


(2.4) 


tefl  +  A  arr),f„i>.  a* 

***?■ 

—  w.-b- 


1=/ 
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(2.5) 


In  the  above  equations  the  operator  is  defined  by 


A 

dt 


4-  ^  $~ 
7 1 


The  prime  denotes  differentiation  with  respect  to  x  and  N  is  the 
maximum  number  of  species.  In  these  equations  the  symbols  have  the 
following  meaning:  u,  v,  and  w  are  the  velocity  components  along  the 
x,  y,  z  coordinates  respectively;  f)  k’/*  ,  cp  are  the  mass  density, 
heat  conductivity,  viscosity  coefficient,  and  specific  heat  coeffi¬ 
cient  at  constant  pressure  of  the  gas  mixture;  (i  =  1,  2.  ...N) 
are  the  binary  diffusion  coefficients.  These  gas  properties  are 
functions  of  the  temperature  T  and  of  the  species  concentration  c^  in 
mass  fraction  units.  The  production  and  loss  of  species  through  dis¬ 
sociation  and  recombination  is  indicated  by  the  term  w£  which  is  the 
mass  rate  of  production  of  the  ith  species  per  unit  volume  for  each 
reaction.  In  general,  the  left  side  of  these  equations  are  the  change 
of  the  quantities  concerned  due  to  convective  processes  and  are  balanced 
by  the  diffusive  effects,  such  as  viscosity,  diffusion  and  heat  con¬ 
duction.  Conservation  of  momentum  normal  to  the  surface  is  satisfied 
by  the  condition  for  the  constancy  of  the  pressure  in  this  direction. 

In  addition  to  these  equations  the  pressure  is  given  by  the  equation  of 
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state  for  the  gas  mixture  and  is  given  by: 


-f>  =.  RpT £ (2.6) 

where  R  is  the  universal  gas  constant,  is  the  species  molecular 
weight.  The  species  of  interest  in  this  study  those  from  dissociated 
and  ionized  air  which  may  numbered  to  10  for  the  neutrals  and  their 
positive  ion. 

The  species  equation  satisfy  a  global  conservation  of  mass  given 
by 

2L  C/ov/i)  + 

27C  j  /12H 


=  o 


(2.7) 


These  equations  together  with  chemical  rate  constants,  thermodynamics, 
and  transport  properties  described  in  Section  IV  are  sufficient  to  de¬ 
termine  the  (4  +  N)  unknowns.  Due  to  the  constancy  of  the  pressure 
through  the  boundary  layer  the  pressure  is  prescribed  by  the  inviscid 
flow  outside  of  the  boundary  layer. 

From  the  problem  as  stated  above  the  flow  at  angle-of-attack  over 
a  rotationally  symmetric  body  such  as  a  sphere-cone  combination  has 
symmetry  with  respect  to  the  plane  containing  the  free  stream  velocity 
vector  and  a  body  meridian  curve.  The  meridian  curve  containing  the 
stagnation  point  of  the  flow  is  a  line  of  symmetry  and  the  velocity 
component  w  normal  to  this  curve  along  the  body  is  zero.  One  notes 
also  that  the  equations  for  a  set  of  wind  oriented  coordinate  con¬ 
structed  at  the  stagnation  point  when  it  is  on  the  sphere  reduces 
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immediately  to  the  same  form  as  that  of  the  zero  angle-of-attack  case. 

With  consideration  of  a  boundary  layer  concept  the  gradients  normal 
to  the  surface  dominate  the  diffusive  processes  due  to  viscosity, 
heat  conduction,  and  species  diffusion.  The  implication  of  this  is 
that  the  differential  system  remains  parabolic  and  the  characteristics 
at  any  point  of  the  flow  are  changed  only  by  upstream  conditions.  How¬ 
ever,  a  complication  of  the  three  dimensionality  appears  in  that  the 
streamlines  do  not  remain  in  the  same  meridian  plane  and  imposes  a 
condition  on  the  computational  aspect  of  the  problem. 

For  a  sphere  cone  body  at  non-zero  angle-of-attack  to  the  free 
stream  the  flow  over  most  of  the  spherical  part  of  the  body  is  symmetric 
about  an  axis  parallel  to  the  free  stream.'  That  is,  the  stagnation 
point  of  the  flow  remains  on  the  sphere  if  the  angle-of-attack  is  less 
than  the  complement  of  the  cone  angle.  This  condition  is  satisfied  for 
the  values  of  angle-of-attack  for  which  the  three  dimensional  laminar 
boundary  layer  equations  are  valid.  The  origin  of  the  coordinates  x, 
y  and  z  is  located  at  the  nose  of  the  sphere  cone  body.  This  is  the 
body  oriented  set  of  coordinates.  At  the  stagnation  point  another  set 
of  coordinates  can  be  constructed  which  is  wind  oriented.  The  spherically 
symmetric  solutions  obtained  for  the  zero  angle-of-attack  case  are  valid 
in  this  set  of  coordinates.  By  this  mean  the  flow  characteristics  for 
the  spherical  portion  of  the  flow  can  be  obtained  by  the  zero  angle-of- 
attack  assumption  and  utilized  as  initial  conditions  for  the  angle-of- 
attack  case. 

A  further  transformation  is  made  to  scale  the  growth  of  the  boundary 
layer  due  to  compressibility  effect.  This  is  given  by 
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(2.8) 


h  _ 


5  =  5  (*  )  >  7=2 

In  this  set  of  coordinates  ^  ^  ^  are  the  transformed  coordinates 
along  and  normal  to  the  meridian  and  ^  is  the  meridional  angle. 
The  application  of  this  transformation  (7)  to  equation  (1)  to  (5) 
lead  to  the  following: 


^  f  ZI'L  if  4.  Vlf  4  ^  ^5  f 
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(2,.?) 


(2.11) 


and  the  species  equation  for  each  species 


oC  f  £  .f.  i*£,'  _  JL  f  ji  ?ct 
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The  function  V  is  given  by  the  global  equation  of  mass  (6) 


(  lO^h) _ /_ 

^4  A  j&i/€ 


oC  WL  2.  ^5av\£  ) 
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Equation  (13)  is  obtained  from  the  conservation  of  mass  (6)  by  noting 
that  the  latter  is  satisfied  by  a  vector  potential  whose  derivatives 
lead  directly  to  the  physical  velocity  components.  A  transformation 
of  the  coordinate  then  leads  to  equation  (13).  In  these  equations  the 
notation  is  as  follows: 


/'=u/ue  ,  -£=  C  fA  >  /  Cf./O  >  = 


f  D.  c,  is  the 


Lewis  number , 


=  /*-  O’  ^ 


T/T  .  The  subscript  e  denotes  values  at 


is  the  Prandtl  number,  ^  . 
the  boundary  layer  edge  and  are  obtained  from  the  inviscid  flow  as 
given  in  Section  III.  In  addition,  the  \/  function  is  a  combination  of 
the  velocity  components  of  the  form: 


(2.14 


In  these  equations  o L  and  ^  remain  arbitrary  function  of  their 
respective  arguments  and  are  assigned. 

2. 2  Method  of  Solution  by  Finite  Differences 

The  method  utilized  for  the  solution  of  the  three  dimensional 
boundary  layer  equations  as  set  forth  in  the  previous  section  is  the 
finite  difference  one.  The  coupled,  nonlinear,  partial  differential 
equations  are  quasi-linear  with  second  order  derivatives  in  one  direction 
only,  namely,  the  direction  y  normal  to  the  surface  of  the  body.  Thus 
these  equations  remain  parabolic  with  non-zero  angle-of-at tack  as  in 
the  two  dimensional  case.  The  three  dimensionality  of  the  flow  at 
angle-of-attack  leads  to  a  cross-flow  velocity  component  which  tends  to 
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sweep  the  flow  around  the  body.  Due  to  this  factor  the  domain  of  dependence 
of  a  given  point  depends  on  the  spanwise  mathematical  characteristics  of  the 
equations  in  a  complicated  way.  One  set  of  characteristics  of  this  set  of 
equations  is  the  streamlines  or  traces  of  the  streamlines  in  planes  parallel 
to  the  surface  of  the  body.  These  lines  have  a  cross-flow  direction  to  the 
main  free  stream  and  carries  disturbances  from  one  part  of  the  flow  to  another. 
Disturbances  are  also  propagated  downstream  by  first  diffusing  normal  from 
the  wall  to  the  boundary  layer  edge  where  they  are  then  carried  downstream  by 
the  external  flow.  The  characteristics  of  the  convective  operator  including 
the  cross-flow  component  lead  to  the  Courant-Fr iedr ich-Lewy  condition  of 

(2.15) 

d*  1A 

This  condition  implies  that  domain  of  the  difference  scheme  contains  that  of 

(12) 

the  differential  system  .  The  implication  of  this  stability  condition  is 

that  difference  schemes  which  are  inherently  stable  such  as  implicit  type  are 

now  conditionally  stable  for  three  dimensional  flow;  the  condition  given  by 

equation  (15)  must  be  satisfied  at  each  point. 

The  difference  scheme  of  Crank-Nic olson  or  the  fully  implicit  one  has 

been  used  in  the  zero  angle-of-attack  case.  There  have  been  a  number  of 

generalization  of  these  schemes  (implicit  and  explicit)  to  consider  the  non- 

(13,14,15,16) 

zero  angle-of-attack  flows  .  The  principal  problem  which  pre¬ 

sents  itself  in  three  dimensional  flow  is  the  conditional  stability  introduced 

(13) 

by  the  angular  momentum  equation.  Krause  has  introduced  a  scheme  where¬ 

by  reverse  cross-flow  can  be  considered.  This  can  be  useful  for  such  cases. 

(14) 

A  difference  scheme  related  to  the  Crank-Nicolson  one  for  two  dimensional 
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(18) 


flow  has  been  utilized  in  this  study.  A  similar  scheme  was  utilized  by  Lew 
for  the  study  of  chemically  reacting  gases  and  solids  with  time  dependency. 

The  scheme  averages  the  derivatives  in  the  normal  direction  and  angular  direc¬ 
tion  in  such  a  way  that  the  algebraic  equations  are  tri-diagonal.  Consider  a 
three  dimensional  rectangular  grid  with  grid  lines  parallelling  the  x,  z,  y 
axes  with  the  indices  1,  m,  n  along  the  x,  z,  and  y  directions  respectively. 
The  represnetative  second  order  equation  for  either  the  velocity,  temperature, 
or  species  concentration  can  be  given  by 

=  F  ($, 5,1 ',£>■$,>%, $s  )  (2-16) 

where  <j)  is  any  of  the  dependent  variables.  The  finite  difference  re¬ 
presentations  are  then 
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These  equations  for  the  difference  representation  lead  to  difference 

equations  involving  the  three  values  of  n  along  two  surface  normals  diagonal 

to  each  other  and  two  known  additional  points.  The  finite  difference 

scheme  leads  to  a  tri-diagonal  scheme  for  three  unknowns  for  each  equation 

at  the  points  along  the  normal  passing  through^  +  1  and  m  +  1.  The 

central  differences  are  used  for  the  first  derivatives.  Any  instability 

in  the  implicit  scheme  arises  from  the  cross- flow  in  the  m  index  direction. 

The  first  derivatives  with  respect  to  the  normal  distance  are  evaluated 

at  diagonal  points  instead  of  at  the  four  vertices  of  the  rectangle  as  pres- 

(13) 

cribed  by  the  Crank-Nicolson  scheme.  Krause  utilizes  a  forward  point 
difference  of  the  cross  derivative  and  can  account  for  reverse  cross-flow. 
This  would  be  a  simple  addition  to  the  present  scheme.  The  tri-diagonal 
algebraic  equation  can  be  written  in  the  form: 

An.  Wn+i  ■+  E>„  Wn  -f  C-n  ^n-i  •=?  (2.17) 

(12) 

and  can  be  solved  by  a  simple  algorithm  .  The  coefficients  An,  Bn  ,  and 
Dn  are  obtained  from  one  sweep  progressing  in  n;  then  the  unknown  wn  is 
obtained  from  a  second  sweep  starting  from  the  maximum  value  of  n.  Thus, 

wn  ^  E„  v\znf/  +  e*  (2.18) 

where 

=—  CAa.+Cs,F  )/ (b^+C^H) 

£„  =  -Ajah  +G,  £„_, ) 

ch  -  CT)n )/ CBn  +C„  £„.,) 


-  16  - 


for  3  n  ^  N-l  where  N  is  the  maximum  number  of  points  used. 

The  boundary  conditions  are  given  in  the  same  form  and  can  be  written  as 

w,  -  +  PvVj  +  -h 

in  terms  of  the  boundary  values  coefficients  of  H,  F,  h,  and  g.  The 
solution  involved  in  the  above  is  valid  for  a  set  of  linear  algebraic 
equations.  The  set  of  nonlinear  differential  equations  from  which  the 
finite  difference  equations  were  obtained  was  linearized  first  by  the 
concept  of  quasi-linearization.  This  quasi-linearization  is  required 
for  the  chemical  production  terms  in  both  the  species  and  energy  equa¬ 
tions  and  in  the  nonlinear  convective  terms. 

2.3  Initial  Conditions 

For  the  sphere  cone  body  with  a  large  nose  radius  such  as  the 
Trailblazer  body  the  stagnation  point  will  locate  itself  on  the  spheri¬ 
cal  portion  of  the  body  for  a  large  range  of  values  of  the  angle-of- 

attack.  The  flow  over  the  sphere  is  symmetric  with  respect  to  the  wind 
coordinates  when  the  stagnation  point  is  located  there.  The  relative 
coordinate  systems  are  shown  on  Figure  2.1.  The  set  of  equations  as  given 
in  Section  2.1  can  also  be  utilized  for  the  flow  with  respect  to  the  wind 
coordinates;  it  reduces,  of  course,  to  the  two  dimensional  flow  (at  zero 
angle-of-attack)  considered  previously  in  reference  (3).  The  wind  coordinate 
symmetric  flow  solution  will  be  utilized  as  initial  conditions  for  the  flow 
at  non-zero  angle-of-attack  at  a  point  where  both  solutions  are  valid.  This 
usually  occurs  at  the  shoulder  juncture  on  the  windward  meridian.  The  inter¬ 
section  of  the  streamlines  from  the  stagnation  point  with  each  meridian  pro¬ 
vides  the  angular  displacement  of  the  velocity  components  u  and  w.  Since  there 
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is  no  variation  of  characteristics  across  streamlines  in  the  symmetric  por¬ 


tion  of  the  flow  the  conditions  at  the  initial  line  are  obtained  by  an  inter¬ 
polation  for  each  point  of  one  set  of  zero  angle-of-attack  solutions.  The 
velocity  in  the  angular  direction  for  the  body  oriented  coordinates  involve 
a  component  of  axial  velocity  dependent  on  the  angle-of-attack  and  body 
location. 
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XII.  Inviscid  Flow  for  Edge  Boundary  Conditions 


As  indicated  in  Section  II  the  edge  values  of  the  flow  velocities, 
temperature,  and  species  concentration  are  required.  These  values  are 
usually  obtained  from  considerations  of  the  inviscid  equations  for  three 
dimensional  flow.  Since  only  the  edge  values  are  needed  approximations 
by  means  of  utilizing  experimental  fits  of  pressure  distribution  or  its 
Newtonian  approximation  have  been  found  adequate  for  the  zero  angle-of- 
attack  case.  For  a  slightly  blunted  conical  body  streamlines  entering 
the  downstream  boundary  layer  edge  are  not  the  body  streamline  and  an 
experimental  fit  of  the  shock  shape  has  been  found  sufficient  to  obtain 
the  effect  of  entropy  gradient  in  the  inviscid  flow.  This  method  of 
obtaining  the  edge  conditions  has  been  utilized  for  the  symmetric  portion 
of  the  flow  (reference  3). 

3.1  Edge  Conditions  at  Surface  of  Body  at  Angle-of-Attack 

The  edge  conditions  for  the  non-symmetr ic  portion  of  the  body  require 

a  more  complicated  analysis  because  of  the  additional  degree-of-freedom. 

For  the  first  step  the  body  streamline  characteristics  are  obtained  by 
utilizing  the  solution  of  the  three  dimensional  equations  ©f  motion  and 
prescribed  pressure  distributions  over  the  entire  body.  The  solution 
involves  the  set  of  momentum,  energy,  and  species  equations  for  a  chemical 
reacting  finite  rate  chemistry  multicomponent  gas  with  dependency  on  the 
two  surfaces  coordinates  x  and  z.  An  iterative  finite  difference  scheme 
has  been  chosen  for  the  method  of  solution. 

The  equations  of  motion  at  the  surface  of  the  body  for  an  inviscid 
flow  can  be  obtained  from  equations  (2,2)  to  (2.5)  by  setting  the  transport 
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terms  to  zero.  At  the  surface,  y  -  o  and  the  normal  component  v  of  the 
velocity  is  also  zero.  Thus  the  operator  d/dt  becom  simply 


d  ^  ]L 


dt 


/  vve  9 

A.  35 


where  the  subscript  e  denotes  the  inviscid  edge  values  for  the  boundary 
layer.  These  equations  are  then 


clVe  —  Vv/Vi/  , 
dt  At 


_  ±  i_£ 

fe  3X 


(3.1) 


chU  ,  V  At  - L 

^  *  £*■  ** 


(3.2) 


cTt  %  Cf  K  3  X  K  ^  l£t 


*5  ct7 


(3.3) 


4fit  -  X  t'e/jA.  .,,A/  . 

St  £  e  ; 


In  addition,  the  equation  of  state  as  given  by  equation  (2.6)  is  used 
to  relate  the  pressure  to  the  other  thermodynamic  variables.  The  enthalpy 
and  chemical  production  terms  are  the  same  as  described  in  Section  2. 

For  conditions  at  the  body  surface  the  pressure  distribution  is  pre¬ 
scribed  by  a  modified  Newtonian  law  dependent  on  the  flow  deflection  at 
the  particular  point.  Reliable  pressure  value  on  the  body  surface  are 
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usually  obtained  from  a  modified  Newtonian  theory  at  high  deflection 
(15) 

angles.  Der  has  suggested  an  improvement  of  the  surface  pressure 

at  the  smaller  deflection  angle  (  45°)  on  the  basis  of  experimental 

data  and  characteristics  solutions  for  a  sphere-cone  (NASA)  and  this  is 
approximated  by  the  equation  of 

~  * +  *if  [+&-£)]  >  °'5> 

Cp  _  COS  G  ;  ^  <-£  (3.6) 

/VfAx 

where  ©  is  the  complement  of  the  flow  deflection  angle 

and  Cp  is  the  coefficient  of  pressure.  Equations  (3.5)  and  (3.6)  are 

matched  at  6  -  TT/4. 

With  the  pressure  distribution  at  each  point  of  the  body  surface 
given  by  equations  (3.5)  and  (3.6)  the  inviscid  edge  flow  can  be  ob¬ 
tained  from  equations  (3.1)  to  (3.4).  There  are  (3  +  N)  unknowns  of 

ue:,  we  the  two  velocity  components,  Te  the  temperature  and  C*  (i=l,2,..N) 

e 

the  species  concentration  in  mass  fraction.  These  nonlinear  partial 
differential  equations  are  first  order  and  require  a  different  method 
of  solution  from  the  previous  section.  This  is  discussed  in  the  next 
section.  The  values  of  the  dependent  variables  are  prescribed  at  the 
initial  line.  This  is  sufficient  to  provide  the  solution  at  any  other 
point.  For  the  study  here  the  initial  line  is  chosen  normal  to  its 
body  symmetry  axis,  that  is,  in  the  body  oriented  coordinate.  As 
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described  in  the  previous  section  this  initial  line  intersects  the 
streamlines  issuing  from  the  stagnation  point  at  different  points 
along  their  flow  (See  Fig.  2.1).  The  velocity  components  at  this 
initial  line  are  given  .in  terms  of  the  wind  oriented  coordinate 
solutions. 

The  relation  between  a  wind-oriented  coordinates  which  has  its 
origin  at  the  stagnation  point  and  a  body-oriented  with  its  origin 
at  front  symmetry  point  can  be  given  by  a  set  of  equations*. 


These  equations  are  valid  for  a  body  at  an  angle-of-attack  to  the 
free  stream  and  x,  s  denote  the  distance  along  the  body  measured  from 
the  body  origin  and  stagnation  point  respectively  and  where  5 ,  t 
denote  the  angular  displacement  with  respect  to  the  corresponding  axes. 

It  is  noted  that  these  quantities  become  equal  wheno6^=  0,  the  zero 
angle-of-attack  case.  Equations  (3.7)  transform  the  initial  profiles 
which  are  symmetric  with  respect  to  the  wind-oriented  axes  into  the  body- 
oriented  coordinates.  The  flow  characteristics  are  symmetric  for  constant 
value  of  s  for  the  wind-oriented  coordinates  and  provide  the  initial  con¬ 
ditions  at  each  body  meridian. 

3.2  Method  of  Solution  of  Inviscid  Edge  Vallues 

The  prescription  of  the  pressure  distribution  by  equations  (3.5)  and 
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(3.6)  with  the  deflection  angle  obtained  from  equations  (3.7)  leads 

directly  to  the  pressure  gradients  required  for  the  solution  of 

equations  (3.1)  to  (3. A).  These  equations  constitute  a  nonlinear 

first  order  set  of  partial  differential  equations  in  two  directions. 

They  are  solved  by  a  predictor-corrector  method  such  as  used  by  Mac— 

(19) 

Cormack  .  In  this  scheme  the  solution  can  be  predicted  by 

F  (m  JL+i )  =  F(yr )  -  A (yn}£)  [ F (n+i )JL) 

(3.8) 

-F(n>*-)]  a$/a$  + 

\ 

and  corrected  by 

F  (  m.,1+1  )  =  -L  I  F  (nijS.  )  +  F(*>1+') 

_/t  fmj  A  h)[  F  (yn1  !+i (3  9) 

This  method  is  useful  for  flows  where  convection  dominates  the  situation 
and  is  stable  provided  the  CFL  condition  is  satisfied.  In  this  case  the 
criterion  for  stability  is  given  by 

A ^  A  J  JA  [TrljJl  )  (3.10) 

and  is  satisfied  by  all  the  equations  at  each  point  of  the  grid.  In 
equations  (3.8)  and  (3.9)  JL  and  m  are  the  indices  along  the  meridian  curve 
measured  by  ^  from  the  stagnation  point  (  1  ■  x  -  Rjj  for  inviscid 
flow  considered  in  this  section)  and  along  the  angular  displacement  measured 
by  5-  The  coefficients  A  and  G  are  given  by  each  individual  equations 
(3.1)  to  (3.4). 
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Initial  conditions  are  required  m  =  1  along  the  stagnation  line 
and  at  a  value  of  JL  £  1.  At  the  stagnation  line  and  the  leeward 
symmetry  line  the  flow  has  zero  angular  velocity  due  to  symmetry. 

Thus  the  velocity  component  we  is  zero  and  equation  (3.1)  to  (3.4) 
reduces  to  flow  along  the  streamline.  These  equations  along  the 
streamlines  are  again  integrated  by  a  predictor-corrector  of  the  same 
form  as  equations  (3.8)  and  (3.9)  except  that  there  is  only  one 
coordinate  involved. 

The  initial  conditions  for  all  the  three  dimensional  equations 
are  obtained  from  the  inviscid  flow  for  the  front  portion  of  the 
sphere 
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IV.  Chemical,  Thermodynamics  and  Transport  Gas  Model 


4. 1  Chemical  Kinetic  Rate  Constants 

The  chemical  kinetics  system  for  a  nitrogen-oxygen  atmosphere  is  of 
concern  for  Trailblazer  II  conditions.  The  flight  regime  involves  prin¬ 
cipally  velocities  17,500  ft/sec.  Since  there  is  no  ablation  of  body 

material  the  gas  mixture  is  composed  of  dissociated  and  ionized  air  spe- 

+ 

ies.  For  most  cases  these  species  are  O2,  N2,  N.  0.  NO,  NO  and  electrons. 

The  higher  temperature  range  to  15,000°K  for  altitudes  around  100  Icm  requires 

+  +  +  + 

the  additional  atomic  and  molecular  ions  of  N2  ,  0  ,  N  ,  and  0  .  The 

chemical  system  involving  the  neutral  species  and  NO  is  given  in  reactions 

1  to  6  in  Table  I  which  also  contains  additional  reactions  which  become 

important  at  the  higher  temperature  range  of  15,000°K.  This  Table  of  re- 

(20) 

actions  has  been  suggested  by  Bortner  for  a  nitrogen-oxygen  system 

typical  of  the  shock  layer  of  re-entry  bodies.  At  a  higher  temperatures 
the  assoicative  ionization  reaction  7  and  8  and  collisional  ionization  re¬ 
actions  for  the  neutral  species  as  given  by  9  to  13  become  important. 

Finally,  after  the  formation  of  the  positive  ions,  exchange  reactions  re¬ 
distribute  these  charges  by  means  of  reactions  14  to  20.  The  reaction 
rate  constants  have  been  suggested  by  Bortner  in  Reference  (21)  and  (22) 
where  he  also  discusses  their  uncertainties  and  limitations.  Results  in 
this  report  are  obtained  using  reactions  1-6  and  reactions  1-10,  14,  15, 

17,  18,  20.  A  minimum  number  of  reactions  consistent  with  reality  is 


utilized . 

In  general,  the  chemical  reactions  are  represented 


by  the  equation 
(4.1) 
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where  the  subscript  r  represents  the  particular  reaction,  i  is  the  species 

that  is  being  reacted  and  cL>  .  lb  .  are  the  stoichiometric  coefficients. 

r  i  '  r 1 

In  this  report  the  values  of  r  and  i  range  over  the  number  of  reactions 
and  the  number  of  species  respectively.  For  this  representation  the  net 
mass  rate  of  production  of  species  i  per  unit  volume  can  be  immediately 
given  by 


fL  f  jl 

W,.=/^;Z.(/3  .-ot  ,)/A  /'  (tS) 

L  r*‘  I  f'l'i  My 


oC  • 


(4.2) 


In  this  equation  is  the  molecular  weight  of  species  1  and  the  subscript 
denotes  the  species  and  the  catalytic  bodies,  i.e,  j  =  1,  2,  ...N,  N+l,  .  .  . 
The  mass  fraction  of  catalytic  bodies  is  given  in  terms  of  the  N  species  by 

J 


=  n-  Z  zr,  u.  . 

3  Cs.i  (3 -A/)/ 


(4.3) 


for  j  =  N  +  1,  ....  J.  The  use  of  equation  (4.3)  facilitates  the  com¬ 
putation  problem  for  large  systems  of  chemical  species  and  flow  variables. 

Equation  (4.2)  is  the  source  term  for  the  mass  generation  of  species 
and  is  a  function  of  the  density  ^  of  the  mixture,  the  temperature  T, 
and  the  density  (  jO  c^)  of  the  pertinent  species.  The  rate  constants 
and  k^r  are,  of  course,  functions  of  the  temperature  T. 

4.2  Thermodynamic 

The  gas  model  is  composed  of  air  species  and  reactions  as  considered 
above.  The  total  enthalpy  of  the  air  mixture  is  easily  obtained  when 
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the  individual  species  values  are  first  made  available.  For  the  chemical 
reacting  gas  mixture,  the  enthalpy  is  obtained  by  the  summation 


Jh  =.  Z  C, 


£  — / 


(4.4) 

(23) 


where  the  pure  species  enthalpy  values  have  been  tabulated  by  Browne 
The  thermodynamic  properties  for  the  diatomic  molecules  were  obtained 
from  the  second  virial  coefficient  approach  which  essentially  leads  to 
the  thermodynamic  properties  departure  from  the  monatomic  gas  at  the 
higher  temperature;  the  Morse  potential  is  utilized.  At  the  lower  temper¬ 
ature,  the  diatomic  species  are  considered  as  rigid  rotator  harmonic  os¬ 
cillators  with  corrections  for  rotational  stretching  and  vibrational  an- 

harmonicity.  Comparison  of  these  data  with  others  are  excellent  and  have 

(23) 


been  given  by  Browne  .  The  individual  species  values  of  the  enthalpy 

are  given  in  Figure  (4.1). 

The  mixture  specific  heat  is  obtained  similarly  as 


(4.5) 


4.3 


Transport  Properties 

The  transport  properties  involve  the  coefficients  of  viscosity  and 

conductivity.  The  individual  species  coefficients  have  been  obtained 

from  a  number  of  sources.  The  species  viscosity  in  terms  of  the  collision 

.  (2,2) 

integrals  is  given  by  the  relation: 


-  6 


CMcT) 

up# 


vz 


fotses 


(4.6) 
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(24) 

Collision  integrals  from  Yun  and  Mason  for  dissociated  air  species 

are  utilized  in  the  compututations ;  the  integrals  are  dependent  on  the 
intermolecular  potential. 

The  conductivity  of  the  species  is  related  to  A  by 


A'=  /HBr  +  M5 
yfc  L  pt  R  j 


where  is  in  g-cal  per  (cm  sec  K). 

Finally,  the  mixture  values  of  the  viscosity  and  conductivity  are  given 
by  Wilke's  rules  whereby 

N 


C  -I 


(4,7) 


\ 


(4.8) 
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Numerical  Results 


Some  numerical  results  are  given  in  this  section;  detailed  flow 

field  calculations  will  be  presented  in  a  separate  data  memorandum. 

These  results  were  obtained  with  the  following  considerations.  The 

six  species  O2 ,  N2,  NO,  N,  0,  and  N0+  with  reactions  1  to  6  were  used. 

For  the  angles-of-attack  of  interest  the  stagnation  point  of  the  flow 

lies  on  the  spherical  portion  of  the  body  resulting  in  a  region  of 

symmetrical  flow.  The  distributions  of  flow  characteristics  for  flows 

with  symmetry  (at  zero  angle-of-attack)  have  been  given  in  the  previous 
(3) 


report  (Lew  ).  These  results  are  used  as  initial  conditions  for  the 
three  dimensional  flow  considered  herein.  One  notes  that  the  three 
dimensional  boundary  layer  program  can  also  be  initiated  directly  at 
the  stagnation  point  by  generating  its  own  intial  profiles  by  appropriate 
use  of  the  proper  coordinates.  By  utilizing  the  two  dimensional  symmet¬ 
rical  flows  as  initial  conditions  there  is  a  slight  saving  in  cost. 
Moreover,  the  effect  of  swallowing  for  the  front  portion  of  the  flow 
is  accounted  for  by  this  means. 

The  pressure  distribution  utilized  in  this  report  is  a  first  approx¬ 
imation.  Previous  results  in  reference  (2)  and  (3)  use  a  next  approx¬ 
imation  to  the  inviscid  edge  conditions.  These  conditions  are  valid  for 
a  thin  shock  layer  with  pressure  variation  across  the  shock  layer  given 
by  the  first  integral  of  the  momentum  equal  normal  to  the  surface.  This 
integral  is  given  by 


■f  pl  M 


(5.1) 


where  pg  (x)  =  p  (x,  "y^  )  Pressure  at  the  shock,  T  is  the 

stream  function{this  may  be  measured  with  respect  to  the  oody  where 
it  has  a  value  of  o)  and  u,  Rc,  r  are  the  tangential  velocity,  radius 
of  curvature,  and  the  local  radius  around  the  axis  of  symmetry  re¬ 
spectively.  Consistent  with  the  approximation  of  this  thin  shock  layer 
theory,  the  integral  is  evaluated  at  the  shock.  This  fact  simplifies 
equation  (5.1)  to  the  form: 

f(P^)^ps  u)-psW)(V-i)  (5.2)  . 

where  the  subscripts  B  and  S  denote  values  evaluated  on  the  body  and 

shock  respectively.  This  relation  was  used  for  the  next  approximation 

(2) 

to  the  pressure  in  the  previous  two  dimensional  results,  see  also 

(25) 

Lew  for  further  discussion  and  derivation.  The  use  of  relation 

(5.2)  require  a  shock  shape  at  each  flight  condition  and  a  number  of 
experimental  and  analytical  results  of  shock  shape  for  hypersonic  flow 
over  sphere-cone  at  zero  angle-of-attack  was  correlated  and  used.  As 
the  results  of  reference  (3)  show,  the  pressure  for  a  blunted  cone  drops 
rapidly  at  the  shoulder  dependent  on  the  cone  angle,  overexpands  below 
the  downstream  conical  values  and  then  rises  asymptotically  to  these  latter 
values.  Dissociated  atomic  species  decay  rapidly  prior  to  the  shoulder 
having  increased  for  a  short  distance  from  the  stagnation  point.  The  species 
concentration  after  the  initial  point  considered  here  are  those  for  the 
body  streamline  after  a  partial  swallowing  has  been  accounted  for. 
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As  examples  the  inviscid  edge  values  of  velocity,  temperature, 
and  species  concentration  were  obtained  fro  the  case  of  the  AFCRL 
Trailblazer  body  at  240,000  ft.  altitude  flying  at  17,500  ft. /sec. 
for  non-zero  values  of  angle-of-attack.  These  values  are  2,  6  and  12 
degrees.  Since  there  is  symmetry  in  the  plane  of  the  velocity  vector 
only  the  variation  from  ^  =  0  to  7T  is  considered.  Initial  condi¬ 
tions  were  obtained  from  reference  (2).  Results  are  given  on  Figures 
(5.1),  (5.2)  and  (5.3)  for  the  cases  of  angles-of-attack  indicated. 

In  Figure  (5.1)  the  variation  of  the  angular  velocity  component  around 
the  body  shows  a  peak  at  about  90  degrees  azimuthal  angle.  This  peak 
increases  with  angle-of-attack.  Thismeans  that  the  streamlines  have 
their  largest  directional  deviation  from  the  axis  of  symmetry  at  this 
point  and  gradually  align  themselves  to  the  leeward  line  of  symmetry  and 

then  proceed  downstream.  The  angular  displacement  of  the  streamline 

-1 

from  the  meridian  plane  is  given  by  tan  (w£/ue).  There  is  a  drop  of 
temperature  in  the  azimuthal  direction  as  expected;  the  pressure  is  lower 
on  the  leeward  side.  The  axial  velocity  changes  little  from  the  windward 
to  leeward  planes.  Figure  (5.2)  shows  the  variation  of  we  with  angle-of- 
attack.  These  results  behave  qualitatively  as  the  sharp  cone  at  non-zero 
angle-of-attack  to  the  flow.  The  atomic  oxygen  distribution  around  the 
body  at  x  =  2.4  ft.  is  given  in  Figure  (5.3)  and  shows  only  a  small  varia¬ 
tion. 
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The  edge  values  of  the  inviscid  flow  have  been  discussed  above. 

These  have  been  utilized  as  boundary  conditions  for  the  three  dimen¬ 
sional  boundary  layer.  The  same  number  of  species  and  reactions  is 
used  in  this  illustrative  example.  Gradients  for  the  edge  characteristics 
from  the  inviscid  flow  are  obtained  by  a  cubic  spline  fit  for  the  stag¬ 
nation  line  and  a  bicubic  spline  fit  for  the  three  dimensional  aspect 
of  the  boundary  layer.  Typical  results  are  shown  in  Figure  (5.4)  and 
(5.5).  The  gradient  of  the  angular  velocity  is  nonzero  along  the  stag¬ 
nation  streamline  and  is  defined  by 


One  notes  that  the  angular  velocity  component  w  is  zero  at  the  stag¬ 
nation  streamline  to  the  first  power  of  the  angular  displacement  jS  . 

Thus  the  ratio  w/we  has  a  nonzero  value  there.  The  quantity  (we/^  ue) 
has  values  which  vary  at  x  =  .48  ft.  from  U.181  for  b°  angle-of- 

attack  to  0.409  for  12°  angle-of-attack.  This  parameter  gives  an  idea 
of  the  increasing  angular  velocity  as  compared  to  the  axial  velocity 
for  increasing  angle-of-attack.  The  variation  of  g*  across  the  boundary 
layer  is  shown  in  Figure  (5.4),  The  slight  variation  for  small  values  of 
x  increases  downstream  to  an  order-of-magnitude  at  x  =  2.55  ft.’  Devia¬ 
tion  of  the  axial  streamlines  from  the  windward  rays  to  the  leeward  side  is 

I 

thus  indicated.  g  is,  of  course,  zero  in  two  dimensional  symmetrical  flow. 
The  electron  density  distribution  for  two  values  of  x  are  shown  on  Figure 
(5.5).  The  effect  of  angle-of-attack  and  the  pressure  distribution  is 
evident  from  the  figure.  It  is  remembered  that  the  boundary  edge  values  for 
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the  three  dimensional  case  are  evaluated  at  the  body  streamline  after  a 
partial  swallowing  has  occurred. 
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VI .  Concluding  Remarks 


The  problem  of  determining  the  electron  density  distribution  in 
the  shock  layer  of  re-entry  bodies  in  hypersonic  flight  for  nonzero 
angle-of-attack  has  been  considered.  This  report  constitutes  Part  I 
of  a  series  of  report  on  three  dimensional  effects.  In  this  first 
part  the  complete  governing  system  of  equations  for  the  flow  of  an 
ionized  and  dissociated  air  three  dimensional  boundary  layer  and  the 
geometrical  considerations  have  been  given  together  with  the  transport 
and  thermodynamics  model,  and  kinetic  rate  constants  utilized.  In 
addition,  equations  are  developed  for  a  given  pressure  distribution 
for  the  determination  of  the  edge  values  required  for  the  solution 
of  the  boundary  layer  equations.  Finite  difference  methods  have  been 
utilized  for  both  the  solution  of  the  boundary  layer  and  the  inviscid 
flow  necessary  for  the  edge  conditions,  These  methods  are  the  implicit 
difference  scheme  and  the  predictor-corrector  scheme  respectively. 

In  addition  to  the  presentation  of  the  model  and  system  of  equations 
this  report  also  contains  results  for  the  three  dimensional  inviscid 
flow  for  several  angles-of-attack  and  the  results  for  the  stagnation 
line  boundary  layer  at  nonzero  angle-of-attack.  Detailed  evaluation 
of  further  results  will  be  given  in  subsequent  reports. 
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VIII.  TABLE 
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TABLE 

CHEMICAL  SYSTEM  OF  AIR  REACTIONS  AND  RATES 
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a  b  c 

BACKWARD 


°2  +MX 

— 

20  +  Mj^ 

8.74 

(17) 

-1 

.594 

(5)  7.26 

(14) 

-  .5 

N2  +M2 

= 

2N  +  M2 

2.41 

(18) 

-  .8 

.1132 

(6)  1.34 

(17) 

-  .8 

NO  +  Mg 

= 

0  +  N  +  M3 

1.45 

(16) 

0 

.754 

(5)  3.63 

(15) 

0 

0  +  NO 

= 

N  +  02 

4.33 

(7) 

1.5 

.181 

(5)  1.81 

(8) 

1.5 

0  +  n2 

= 

N  +  NO 

6.02 

(12) 

0 

.38 

(5)  1.32 

(13) 

0 

0  +  N 

= 

N0+  +  e 

5.18 

| 

(ID 

0 

.319 

(5)  1.44 

(21) 

-1.5 

N  +  N 

= 

N2+  +  e 

4.7 

(8) 

1.3 

.678 

(5)  5.3 

(17) 

-  .2 

0  +  0 

= 

02+  +  e 

6.62 

(7) 

1.3 

.807 

(5)  7.22 

(18) 

-  .7 

N  +  M 

= 

N+  +  e  +  M 

1.81 

(19) 

-1. 

.169 

(6)  1.09 

(27) 

-2.5 

0  +  M 

= 

0+  +  e  +  M 

3.01 

(18) 

-1. 

.153 

(6)  7.25 

(26) 

-2.5 

N2  +  M 

= 

N  +  e  +  M 

2 

1.81 

(4) 

2 

.181 

(6)  1.09 

(12) 

.  5 

02  +  M 

= 

02+  +  e  +  M 

6.02 

(3) 

1.5 

.14 

(6)  5.8 

(ID 

0 

NO  +  M 

= 

N0+  +  e  +  M 

6.02 

(17) 

-1. 

. 1079' 

(6)  3.63 

(26) 

-2.5 

»2+  -  02 

,  = 

N  +  0  + 

2  2 

1.20 

(14) 

0 

0 

3.61 

03) 

0 

n9+  +  0 

= 

N0+  -+  N 

1.51 

(14) 

0 

0 

£.71 

04) 

0 

0 

0 

0 

.3  (4) 

0 
0 
0 
0 
0 
0 
0 
0 
0 

.409  | 5) 
.354  (5) 


TABLE 


CHEMICAL  SYSTEM  OF  AIR  REACTIONS  AND  RATES  (Continued) 
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FORWARD 
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CHEMICAL  SYSTEM  OF  AIR  REACTIONS  AND  RATES  (Continued) 


Third  Body  Efficiencies 
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FIGURE  4.1  SPECIES  ENTHALPY 
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FIGURE  5.1  AZIMUTHAL  DISTRIBUTION  OF  INVISCID  VELOCITY  COMPONENTS,  TEMPERATURE  AND  PRESSURE 
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FIGURE  5.3  AZIMUTHAL  DISTRIBUTION  OF  ATOMIC  OXYGEN  AT  THE  BOUNDARY  LAYER  EDGE 
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FIGURE  5.5  ELECTRON  DENSITY  DISTRIBUTION  AT  THE  WINDWARD  STAGNATION 
LINE  FOR  ANGLE -OF-ATTACK 


